Laser-assisted electron-argon scattering at small angles 
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Electron-argon scattering in the presence of a linearly polarized, low frequency laser field is studied 
theoretically. The scattering geometries of interest are small angles where momentum transfer is 
nearly perpendicular to the field, which is where the Kroll- Watson approximation has the potential 
to break down. The Floquet R matrix method solves the velocity gauge Schrodinger equation, using 
a larger reaction volume than previous treatments in order to carefully assess the importance of the 
long range polarization potential to the cross section. A comparison of the cross sections calculated 
with the target potential fully included inside 20 and 100 a.u. shows no appreciable differences, which 
demonstrates that the long range interaction can not account for the high cross sections measured 
in experiments. 



I. INTRODUCTION 

Since the development of intense lasers, investigated 
the modification of familiar processes by the presence 
of coherent light. Particles that collide in the presence 
of a laser field can also exchange energy with the field 
in multiples of the laser frequency, a process known as 
laser-assisted collision. This paper investigates free-free 
transitions in electron scattering, where the state of the 
target atom remains unaltered. 

In 1973 Kroll and Watson derived an expression for 
the cross section of a laser assisted scattering event in 
terms of the elastic field-free cross section. Their main 
assumptions were that the time duration of the interac- 
tion is short compared to the laser cycle, and that the 
target itself is unperturbed by the laser. They found the 
following expression, [l| 
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where v denotes the number of laser photons absorbed by 
the electron (— v is the number of photons emitted), kf 
and ki are the momenta of the final and initial electron, 
J„ is the Bessel function of the first kind of order v, 
and is the field free elastic cross section. With the 
time dependent laser vector potential of the form A(i) = 
Aq sin(wt), the other parameters are defined as follows: 
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In addition to assuming that the laser frequency uj is 
small, the Kroll- Watson approximation (KWA) assumes 
that the dimensionless quantity 



mcoj 
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is sufficiently large that only on-shell scattering con- 
tributesQ This assumption becomes questionable not 
just when the parameters of the laser are changed, 
but also in certain critical scattering geometries where 
A • Q « 0. 

The group of Wallbank and Holmes has performed 
experiments investigating these geometries for several 
neutral targets, beginning with helium and argon 0- 
4]. Their results show cross sections orders of magni- 
tude greater than the Kroll- Watson prediction in regions 
where — — Aq • Q is small. One of the first ideas pro- 
posed to explain this discrepancy was the polarization of 
the target by the field of the laser, however several sep- 
arate theoretical treatments 0-0] showed such an effect 
to be negligible. It was also shown that for certain den- 
sities of the target gas, double scattering could account 
for the experimental signal Q. Their determination of 
the experimental density was approximate, however, and 
to our knowledge the density dependence has not been 
confirmed experimentally. 

Madsen and Taulbjerg [9| explore the derivation of the 
KWA, and they develop a generalized approximation by 
expanding the T matrix for weak fields and soft photons, 
but without assuming ([3]) is large. The region where the 
KWA loses its validity is therefore avoided. Their calcu- 
lations show a few scattering geometries where the dif- 
ferential cross section is comparable to the experiments. 
The shape of the experimental cross sections disagrees 
with the Ref. [|| theory, though, and there would have to 
be a very large uncertainty in the determination of the 
elecron scattering angle for their calculations to explain 
the experimental cross sections at all angles. 

The above treatments assume, as does the KWA, that 
only on-shell terms, i.e. terms where energy is conserved, 
contribute to the scattering event. A few later works use 
approximations that include off-shell contributions. Sun 
et al. [10] apply the second Born approximation. Jaron 
and Kaminski jTl], [l2[ also use a similar off-shell approx- 
imation. They suggest that a diffraction effect due to a 
long range potential, i.e. an interaction with extent large 
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compared to the electron deBroglie wavelength, could 
give rise to the sorts of cross sections at small angles 
seen in the experiment. The results of these in compari- 
son with our calculations are discussed further in section 

El 

The advantage of R matrix methods is that they pro- 
vide an exact solution of the Schrodinger equation within 
the chosen reaction volume, and so are limited only by 
the size of that volume and the numerical methods used 
in the calculation. Chen and Robicheaux (l3j use a mixed 
gauge R matrix method with a reaction volume of 30 a.u. 
They calculated cross sections with similar order of mag- 
nitude to the KWA. 

The goal of this work is to use an exact R matrix 
method with an expanded reaction volume. Including 
a longer range for the electron to interact with the in- 
duced dipole potential of the argon atom will clarify what 
contribution the long range interaction has to the laser 
assisted cross section. In section [TT] the details of the Flo- 
quet R matrix method in the velocity gauge are laid out, 
and in section HTT1 the connection to scattering states and 
the form of the scattering matrix are derived. Section HVl 
discusses our numerical results, and section [V] summa- 
rizes our findings. Atomic units are used throughout the 
rest of this paper. 



II. FLOQUET R MATRIX METHOD 

The time and angle dependence of the electron wave- 
function is represented by expanding in a product set of 
spherical harmonic and Floquet basis functions. 
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Here /3 enumerates a complete set of linearly indepen- 
dent solutions of the Schrodinger equation. In order to 
make the treatment of larger reaction volumes tractable, 
this paper treats only scattering geometries where the in- 
coming electron is parallel to the laser polarization. The 
cylindrical symmetry of this case allows setting m z = 
for the entire calculation, considerably reducing the num- 
ber of basis functions needed. 

Taking (3 as a column index, and combining the others 
into a row index, the radial functions can be thought of 
as comprising a matrix F_(r). With this the logarithmic 
derivative, or R matrix, is defined as 



R(r) = F{r) (F'(r))~ 



(7) 



The R matrix describes the behavior of the channel func- 
tions at the surface of a volume of constant radius r. 

The R matrix is found by solving the Hamiltonian in 
the velocity gauge for an electron scattering off of a po- 
tential. The wavelength of the laser is around 2 x 10 5 
a.u., which is much larger than the region of interaction 
of interest, so the vector potential is essentially constant 



in space. This leads to the following form for the Hamil- 
tonian: 



H = \ - -A • p + V(r) 
I c 
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The target atom is represented by a model potential, 
borrowed from Chen and Robicheaux, [l3l j containing a 
shielded Coulombic core and a long range induced-dipole 
term. 
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Here Z = 18, the atomic number, and a — 10.77, the 
argon polarizability. The other parameters of the model 
potential, which were fitted to the field free phaseshifts 
for argon, are ai = 3.04, a 2 = 10.62, 0,3 = 1.83, r cut = 
1.76. 



A. Variational principle for the R matrix 

An extension of the eigenchannel R matrix method, 
adapted for the Floquet formalism and for the veloc- 
ity gauge, yields the solution to the Schrodinger equa- 
tion. The solution is calculated numerically within a fi- 
nite reaction volume S. It is helpful to begin with the 
Schrodinger equation for the velocity gauge hamiltonian 
in integral form. (Note that, for notational brevity, we 
employ notation commonly used in differential geometry 
throughout this section. Function arguments and differ- 
entials are omitted from the integrands, but the integral 
is unambiguous as the domain of integration is denoted 
as a subscript on the integral sign.) 
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The full version of the energy operator is necessary due 
to the fact that, in the Floquet formalism, wavefunctions 
can be superpositions of different pseudoenergy states. 
The only restriction on the space of wavefunctions con- 
sidered is that the spatial inner product between any two 
wavefunctions is periodic, i.e. 



(11) 



This is a reasonable assumption based on the fact that 
the Hamiltonian is itself periodic. 

Application of the first Green identity to the kinetic 
term gives a term with the derivative on the surface <9£. 
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This is the usual starting point to find a variational prin- 
ciple for the logarithmic derivative. Because the velocity 
gauge contains a first derivative term, however, this must 
be taken into account in order to construct a variational 
principle. Using a form of the divergence theorem, the 
identity becomes 



1 
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Now define the operator L as follows: 



— + -A • n* 
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The set of channel functions in all surface coordinates 
forms a linear space on the surface of S. Consider wave- 
functions p that are eigenfunctions of L on this surface, 



i.e. L 9p\ gi , = bp ^| as • The value bp is the usual log- 
arithmic derivative with one term added. 
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The following functional is the variational principle for 
the eigenvalues bp of the generalized logarithmic deriva- 
tive operator: 



bm = 2- 



(16) 

It can be shown that the first variation (^[^J vanishes 
for all deviations ^VP from the exact solution. To show 
this, integrate by parts and use the periodicity restriction 
(ITT1) on the energy operator. 
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Using this, the variation can be written 
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The first Green identity can now be applied to each of 
the kinetic terms, and the divergence theorem can be ap- 
plied to the vector potential terms, to show the variation 
vanishes. 



B. Solving for the R matrix 

The numerical solution of the R matrix is carried out 
by expanding the wavefunction in a basis set: ^> = 

J2 pi c pi^pi( r ^) = J2pi Cpi^T^^ii^it), The radial basis 
functions u p (r) can in principle be arbitrary. In this cal- 
culation we have chosen a finite clement DVR basis, of 
the kind described in [l4j]. The channel functions <I?i(f2, t) 
have the form $j = YuoOCl)- — , + * * . The variational 

principle (|16[) is then written as an eigenvalue equation, 
Ac6 = Tc . 

T pi , qj = r pi (h - i^j ^ + J gsT r p Mv (19) 
Vaj = / V&^ai (20) 

JdS,T 

Here L is the usual Bloch operator with the field term, 
defined so that L^S = + 1A • nVff whose eigenvalues 

r or c 5 ° 



are the logarithmic derivative of the reduced wavefunc- 
tion. This differs from L in equation (1141) only by the 
addition of one hermitian term, so it is still variational. 
Only a small subset of the basis functions have nonzero 
value on the surface of the reaction volume. Denoting 
these open type basis functions with o and the others 
with c, the eigenvalue equation can be written as follows. 

The eigenvalue equation can then be rearranged as fol- 
lows. 

a 00 c 6 = £c n = r 00 -r oc r c " c 1 r co (22) 

This places most of the load of the calculation on the 
linear solution for r^. r co , which requires fewer resources 
than a full diagonalization, although it must be solved for 
each collision energy of interest. The matrix R is then 
calculated from the eigenvalues and eigenvectors found 
in the above. The resulting R matrix is symmetric, and 
is related to the reduced wavefunction F_(r) as follows: 

R-\r) = F'{r) (F(r)y 1 + W (23) 

Wij = - f $?A ■ (24) 
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III. SCATTERING IN A LASER FIELD 
A. Matching to scattering solutions 

The vector potential term in the Schrodinger equation 
can be removed by the following transformation [15j : 



$ = £1* il = exp (- J A(t) • p dr 



(25) 



After defining a{t) = - J A(r) dr, it becomes clear that 
£1 is a translation operator, i.e. for functions of position, 
fi/(r) = /(r + a(t))Q,. It follows simply that $ obeys 
the equation 



^--^V 2 $ + ^(|r + a(t)|) 



(26) 



which at large r approaches the free space Hamiltonian, 
because V(r) — > faster than -, We may therefore match 
to free space scattering solutions, of the form 



-i(E+Viu)t 



(27) 



Here K_ is the short range reaction matrix and the scat- 
tering states are 



fi = 



Iky 



■rn h (k Uz r) 
(28) 

Using the reverse translation, ^ = Q x $, and defining 
functions p(t) and 6{t) that describe the length and angle 
of p(t) = r — a{t) the velocity gauge wavefunction can 
be written as 



P(t) 



(29) 



Following Varro and Ehlotzky [T(|, this can be projected 
onto the original basis set in the untranslated coordi- 
nates, resulting in 



tig 
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(30) 
(31) 
(32) 
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The wavefunction can then be expressed in matrix form 
for radial functions, 



M(r) = f(r) - Q(r)K 



(34) 



and using (|23)) . K can be found in terms of R. 

K = {{R- 1 (r Q )~W)Q(r ) -fl^o))" 1 

((R-\r )-W)l(r ) -R(r )f(r )) (35) 

Note that although W_, f, and g are not symmetric, this 
approach yields a reaction matrix that is symmetric and 
real. 

B. Calculating the cross section 

The asymptotic form of the wavefunction in this Flo- 
quet picture is 



^/ (r,t) = e l(kttZ - Et) + 



j,{k„r-(E+vu))t) 



(36) 



The cross section for each Floquet channel follows from 
this. 

da v | r • jout,iv | fo, . . ..2 , , 

~ln = il — i = T - \J V \")\ (3') 

dil |jincj k 

This can be expressed in terms of the scattering matrix 
S = {1 + iK) {I- iK)~ x . 

do u 
~dtt ~ 
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Y, V* (2i' + 1) (i l '- l S v i,oi> ~ 5i,vS vfi ) Y lfi {6, 0) 



Ll'=0 



(38) 



The phase factor i l 1 above is a result of the choice of 
phase in the scattering functions 



IV. RESULTS 

Figure []] contains our calculation of the cross sec- 
tion along with a comparison calculation using the Kroll 
Watson approximation. The collision parameters were 
chosen to mimic the experiment: a laser intensity of 
5 x 10 7 W cm~ 2 , photon energy of 0.12 eV, and elec- 
tron energy 10 eV. It is worth noting that the electron 
energy lies below the first excitation channel of argon, 
so it makes sense to keep just the single channel for the 
target. The boundary of the reaction volume for this cal- 
culation is at ro = 100 a.u., and it contains 19 Floquet 
channels up to v = ±9 and angular momentum channels 
up to I — 150 for a total of 2869 channels on the sur- 
face of the volume. This number of angular momenta is 
needed to converge the matching equations for the scat- 
tering solution; as is shown in equation (|3ip . the regular 
scattering solutions are proportional to ji (kro), which has 
significant value when I kro- The number of Floquet 
channels needed for convergence increases with rg as well, 
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FIG. 1. (color online) Differential cross section for electrons 
absorbing (top) to 2 (bottom) photons. The solid red line is 
the cross section found with the Floquet R matrix calculation, 
and the dashed blue line is calulated using the Kroll Watson 
approximation. This calculation was performed with the R 
matrix boundary at 100 a.u., with 19 total Floquet channels 
and angular momentum channels up to I — 150. 



but this has more to do with the numerical convergence of 
the R matrix solution than with the matching. The field 
is constant in space, so a larger reaction volume means 
more matrix elements coupling the Floquet channels. For 
Floquet channels v = 1 and v — 2, the differential cross 
sections are converged with respect to basis size to within 
10~ 2 a.u. and 10 -3 a.u. in absolute units of the maxi- 
mum difference between calculations, and to within 1% 
of the cross section value at all angles. 
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FIG. 2. (color online) Differential cross sections at 9 degrees 
versus the number v of photon energies gained by the electron. 
The cross section found via the R matrix calculation (blue 
dots) has comparable order of magnitude to the Kroll Watson 
approximation (red crosses) for a few photon numbers, and it 
is several orders of magnitude smaller than the experimental 
results for all nonzero v. The parameters of the calculation 
are the same as for figure ffj 



Note that the Floquet R matrix calculation and the 
KWA agree very well for all but very small scattering 
angles. The zero photon cross section is also quite close 
to the experimental elastic electron argon cross section 
found by Furst et al. pj}, as expected. Figure [5] shows 
the differential cross section at 9 degrees for absorption 
and emission of up to 4 photons, again compared with 
the KWA result. Note that though they do not agree 
exactly, they differ with far fewer orders of magnitude 
than measured by Wallbank and Holmes at this scatter- 
ing angle. For example, the differential cross sections 
they found for exchanges of 1 and 2 photons were on the 
order of 1 percent of the field free clastic differential cross 
section, while our calulations show these as roughly 10 -2 
percent and 10 -6 percent respectively. 

Differential cross sections calculated with R matrix 
boundaries from 10 a.u. to 100 a.u. show no differences 
that are distinguishable from the convergence with re- 
spect to the basis. Our differential cross section also 
agrees with that of Chen and Robicheaux T3J, who used 
a variable gauge approach and chose an R matrix bound- 
ary of 30 a.u. We can, therefore, rule out the possibility 
that the long range induced-dipole potential would yield 
the sort of diffraction suggested by Jaroh and Kamihski 
jT3 over distances comparable to or even several times 
the electron de Broglie wavelength. 

It has been suggested that uncertainty in the scatter- 
ing angle could account for higher observed cross sec- 
tions. Madsen [9[ even suggests that the incoming elec- 
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tron beam is poorly collimated, leading to an effective 
error in the scattering angle as high as 8°, as opposed to 
the 2° reported from the detector width [4| . Convolving 
the cross sections shown in figure [T] with a gaussian hav- 
ing a width up to 8° does not give a significant difference 
in the 1- and 2-photon cross sections, however. Whatever 
the source of such an error, uncertanties in the scattering 
angle would not explain the experimental cross sections 
for this geometry. 

Sun et ah pi| calculate a cross section for one photon 
exchange that is quantitatively quite close to our result. 
Their result for two photon exchange is several orders 
of magnitude higher, however. This may be due to a 
convergence issue, as a group using a similar method for 
laser-assisted helium scattering at first found high cross 
sections pjj , but later found better results that are closer 
to the KWA pj|. 

V. CONCLUSIONS 

The Floquet R matrix method provides an exact so- 
lution of the Scrodinger equation in the velocity gauge. 



By comparing cross sections calculated with R matrix 
boundaries up to 100 a.u., over ten times the de Broglie 
wavelength of a 10 eV electron, we have shown that the 
induced dipole potential for argon does not contribute to 
the laser assisted cross section. This is true even at small 
angles where the momentum transfer has a very small 
component along the field, and so the KWA is less valid. 
Diffraction from this long range piece of the potential 
can not account for the high cross sections found in the 
experiments of Wallbank and Holmes, which are several 
orders of magnitude higher than both the approximation 
and our calculations. 

The most plausible explanation for the experimental 
results remains multiple scattering. Later experiments 
for helium by the same group claim to see the same high 
cross sections even when the gas is too dilute for multiple 
scattering to play a significant role pl|, so it is unclear 
whether this is the correct explanation. To our knowl- 
edge, a similar experiment for argon including character- 
ization of the gas density has not been performed. 

Discussions with A. Jaron-Becker and M. Tarana are 
much appreciated. We gratefully acknowledge the sup- 
port of the Department of Energy, Office of Science. 
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